Fast and accessible T2 mapping using off-resonance corrected DESPOT2 with application to 3D prostate.

PURPOSE: Most T1 and T2 mapping take long acquisitions or needs specialized sequences not widely accessible on clinical scanners. An available solution is DESPOT1/T2 (Driven equilibrium single pulse observation of T1/T2). DESPOT1/T2 uses Spoiled gradient-echo (SPGR) and balanced Steady-State Free Precession (bSSFP) sequences, offering an accessible and reliable way for 3D accelerated T1/T2 mapping. However, bSSFP is prone to off-resonance artifacts, limiting the application of DESPOT2 in regions with high susceptibility contrasts, like the prostate. Our proposal, DESPO+, employs the full bSSFP and SPGR models with a dictionary-based method to reconstruct 3D T1/T2 maps in the prostate region without off-resonance banding. METHODS: DESPO+ modifies the bSSFP acquisition of the original variable flip angle DESPOT2. DESPO+ uses variable repetition and echo times, employing a dictionary-based method of the full bSSFP and SPGR models to reconstruct T1, T2, and Proton Density (PD) simultaneously. The proposed DESPO+ method underwent testing through simulations, T1/T2 phantoms, and on fourteen healthy subjects. RESULTS: The results reveal a significant reduction in T2 map banding artifacts compared to the original DESPOT2 method. DESPO+ approach reduced T2 errors by up to seven times compared to DESPOT2 in simulations and phantom experiments. We also synthesized in-vivo T1-weighted/T2-weighted images from the acquired maps using a spin-echo model to verify the map's quality when lacking a reference. For in-vivo imaging, the synthesized images closely resemble those from the clinical MRI protocol, reducing scan time by around 50% compared to traditional spin-echo T1-weighted/T2-weighted acquisitions. CONCLUSION: DESPO+ provides an off-resonance insensitive and clinically available solution, enabling high-resolution 3D T1/T2 mapping and synthesized T1-weighted/T2-weighted images for the entire prostate, all achieved within a short scan time of 3.6 min, similar to DESPOT1/T2.

KomaMRI.jl: An open-source framework for general MRI simulations with GPU acceleration.

PURPOSE: To develop an open-source, high-performance, easy-to-use, extensible, cross-platform, and general MRI simulation framework (Koma). METHODS: Koma was developed using the Julia programming language. Like other MRI simulators, it solves the Bloch equations with CPU and GPU parallelization. The inputs are the scanner parameters, the phantom, and the pulse sequence that is Pulseq-compatible. The raw data is stored in the ISMRMRD format. For the reconstruction, MRIReco.jl is used. A graphical user interface utilizing web technologies was also designed. Two types of experiments were performed: one to compare the quality of the results and the execution speed, and the second to compare its usability. Finally, the use of Koma in quantitative imaging was demonstrated by simulating Magnetic Resonance Fingerprinting (MRF) acquisitions. RESULTS: Koma was compared to two well-known open-source MRI simulators, JEMRIS and MRiLab. Highly accurate results (with mean absolute differences below 0.1% compared to JEMRIS) and better GPU performance than MRiLab were demonstrated. In an experiment with students, Koma was proved to be easy to use, eight times faster on personal computers than JEMRIS, and 65% of test subjects recommended it. The potential for designing acquisition and reconstruction techniques was also shown through the simulation of MRF acquisitions, with conclusions that agree with the literature. CONCLUSIONS: Koma's speed and flexibility have the potential to make simulations more accessible for education and research. Koma is expected to be used for designing and testing novel pulse sequences before implementing them in the scanner with Pulseq files, and for creating synthetic data to train machine learning models.

Accurate characterization of dynamic microbial gene expression and growth rate profiles.

Genetic circuits are subject to variability due to cellular and compositional contexts. Cells face changing internal states and environments, the cellular context, to which they sense and respond by changing their gene expression and growth rates. Furthermore, each gene in a genetic circuit operates in a compositional context of genes which may interact with each other and the host cell in complex ways. The context of genetic circuits can, therefore, change gene expression and growth rates, and measuring their dynamics is essential to understanding natural and synthetic regulatory networks that give rise to functional phenotypes. However, reconstruction of microbial gene expression and growth rate profiles from typical noisy measurements of cell populations is difficult due to the effects of noise at low cell densities among other factors. We present here a method for the estimation of dynamic microbial gene expression rates and growth rates from noisy measurement data. Compared to the current state-of-the-art, our method significantly reduced the mean squared error of reconstructions from simulated data of growth and gene expression rates, improving the estimation of timing and magnitude of relevant shapes of profiles. We applied our method to characterize a triple-reporter plasmid library combining multiple transcription units in different compositional and cellular contexts in Escherichia coli. Our analysis reveals cellular and compositional context effects on microbial growth and gene expression rate dynamics and suggests a method for the dynamic ratiometric characterization of constitutive promoters relative to an in vivo reference.

DeepSPIO: Super Paramagnetic Iron Oxide Particle Quantification Using Deep Learning in Magnetic Resonance Imaging.

The susceptibility of super paramagnetic iron oxide (SPIO) particles makes them a useful contrast agent for different purposes in MRI. These particles are typically quantified with relaxometry or by measuring the inhomogeneities they produced. These methods rely on the phase, which is unreliable for high concentrations. We present in this study a novel Deep Learning method to quantify the SPIO concentration distribution. We acquired the data with a new sequence called View Line in which the field map information is encoded in the geometry of the image. The novelty of our network is that it uses residual blocks as the bottleneck and multiple decoders to improve the gradient flow in the network. Each decoder predicts a different part of the wavelet decomposition of the concentration map. This decomposition improves the estimation of the concentration, and also it accelerates the convergence of the model. We tested our SPIO concentration reconstruction technique with simulated images and data from actual scans from phantoms. The simulations were done using images from the IXI dataset, and the phantoms consisted of plastic cylinders containing agar with SPIO particles at different concentrations. In both experiments, the model was able to quantify the distribution accurately.

A Spatial Off-Resonance Correction in Spirals for Magnetic Resonance Fingerprinting.

In MR Fingerprinting (MRF), balanced Steady-State Free Precession (bSSFP) has advantages over unbalanced SSFP because it retains the spin history achieving a higher signal-to-noise ratio (SNR) and scan efficiency. However, bSSFP-MRF is not frequently used because it is sensitive to off-resonance, producing artifacts and blurring, and affecting the parametric map quality. Here we propose a novel Spatial Off-resonance Correction (SOC) approach for reducing these artifacts in bSSFP-MRF with spiral trajectories. SOC-MRF uses each pixel's Point Spread Function to create system matrices that encode both off-resonance and gridding effects. We iteratively compute the inverse of these matrices to reduce the artifacts. We evaluated the proposed method using brain simulations and actual MRF acquisitions of a standardized T1/T2 phantom and five healthy subjects. The results show that the off-resonance distortions in T1/T2 maps were considerably reduced using SOC-MRF. For T2, the Normalized Root Mean Square Error (NRMSE) was reduced from 17.3 to 8.3% (simulations) and from 35.1 to 14.9% (phantom). For T1, the NRMS was reduced from 14.7 to 7.7% (simulations) and from 17.7 to 6.7% (phantom). For in-vivo, the mean and standard deviation in different ROI in white and gray matter were significantly improved. For example, SOC-MRF estimated an average T2 for white matter of 77ms (the ground truth was 74ms) versus 50 ms of MRF. For the same example the standard deviation was reduced from 18 ms to 6ms. The corrections achieved with the proposed SOC-MRF may expand the potential applications of bSSFP-MRF, taking advantage of its better SNR property.

MAPL1: q-space reconstruction using ℓ1 -regularized mean apparent propagator.

PURPOSE: To improve the quality of mean apparent propagator (MAP) reconstruction from a limited number of q-space samples. METHODS: We implement an ℓ1 -regularised MAP (MAPL1) to consider higher order basis functions and to improve the fit without increasing the number of q-space samples. We compare MAPL1 with the least-squares optimization subject to non-negativity (MAP), and the Laplacian-regularized MAP (MAPL). We use simulations of crossing fibers and compute the normalized mean squared error (NMSE) and the Pearson's correlation coefficient to evaluate the reconstruction quality in q-space. We also compare coefficient-based diffusion indices in the simulations and in in vivo data. RESULTS: Results indicate that MAPL1 improves NMSE in 1 to 3% when compared to MAP or MAPL in a high undersampling regime. Additionally, MAPL1 produces more reproducible and accurate results for all sampling rates when there are enough basis functions to meet the sparsity criterion for the regularizer. These improved reconstructions also produce better coefficient-based diffusion indices for in vivo data. CONCLUSIONS: Adding an ℓ1 regularizer to MAP allows the use of more basis functions and a better fit without increasing the number of q-space samples. The impact of our research is that a complete diffusion spectrum can be reconstructed from an acquisition time very similar to a diffusion tensor imaging protocol.